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Abstract 

We explore two Bi-large Neutrino Mixing anzatse within the context of Abelian 
flavor symmetry theories: (BLi) sin 0x2 ~ A, sin #13 ~ A, sin #23 ~ A and (BL2) sin #12 ~ 
A, sin ~ A, sin 023 ~ 1 — A. The first pattern was proposed by two of us and is 
favored if the atmospheric mixing angle 023 hes in the first octant, while the second 
one is preferred for the second octant of 023- In order to reproduce the second texture, 
we find that the flavor symmetry should be U(l) x Z m , while for the first pattern one 
the flavor symmetry should be extended to U(l) x Z m x Z n with m and n of different 
parity. Explicit models for both mixing patterns are constructed based on the flavor 
symmetries U(l) x Z3 x Z4 and U(l) x Zi- The models are extended to the quark 
sector within the framework of SU(5) grand unified theory in order to give a successful 
description of quark and lepton masses and mixing simultaneously. Phenomenological 
implications are discussed. 



1 Introduction 



Our knowledge of the neutrino oscillation parameters has enormously improved in recent 
years. In particular the Daya Bay PQ, RENO [2] and Double Chooz [3] collaborations have 
established that the reactor mixing angle 6* 13 > at about 5 confidence level, confirming 
the early hints for a non-zero 6*13 jUEj. Recent global analyses [6l[7] of neutrino oscillation 
parameters, including the data released at the Neutrino-2012 conference, find that #13 is non- 
zero at about 10cr, and non-maximal atmospheric mixing angle #23 is preferred. However, it 
still isn't clear which octant 623 lies in. The global fit of Ref. [6], prefers 823 in the second 
octant with the best fit value sin 2 623 = 0.613 (0.600) for normal (inverted) neutrino mass 
hierarchy, although this hint is quite marginal and first octant values of 623 are well-inside 
the la range for normal hierarchy and at 1.2(7 for inverted spectrum. While the independent 
phenomenological analyses of atmospheric neutrino data in Refs. [7] obtains a preference 
for 623 in the first octant for both mass hierarchies and exclude maximal mixing at the 2a 
level, the best fit value is found to be sin 2 623 = 0.386 (0.392) for normal (inverted) neutrino 
spectrum. Alternative recent global fits claims both the first and second 823 octants are 
possible [8]. As for the mass-squared difference, the best fit values of Am s 2 ol and Am 2 tm 
are 7.62 x 10~ 5 eV 2 and 2.55(2.43) x 10~ 3 eV 2 respectively, which leads to Am 2 ol /Am 2 tm ~ 
0.030(0.031). Here the values shown in parentheses correspond to the inverted neutrino mass 
hierarchy. Note that the three groups give almost the same 3a ranges for the lepton mixing 
parameters. 

From the theoretical or model building point of view, one implication of these significant 
experimental progress is that it excludes the tri-bimaximal mixing ansatz for neutrino mixing 
[9], unless the underlying theory is capable of providing sufficiently large corrections. So far 
many suggestions have been advanced to explain the new data in particular the largish 

013 [?,man[i5]. 

Instead of seeking for new mass-independent lepton mixing matrices to replace the tri- 
bimaximal pattern [TTHT3] . which may be derived from certain discrete flavor symmetries, 
Ref. jl4j proposed a novel Wolfenstein-like ansatz for the neutrino mixing matrix. In this 
scheme, all the three lepton mixing angles are assumed to be of the same order to first 
approximation 

sin6'i2~A, sin #13 ~ A, sin6 l 2 3~A, (1) 

where A ~ 0.23 is the Cabibbo angle, and the symbol "~" implies that the above relations 
contain unknown factors of order one, the freedom in these factors can be used to obtain 
an adequate description of the neutrino mixing. Inspecting the global data fitting [SHE], 
we see that sin 6*12 — 2.5A and sin 6*13 ~ A/a/2, which is proposed in the so-called Tri- 
bimaximal-Cabibbo mixing [15] and also appeared in the in the context of quark-lepton 
complementarity [16] . Such bi-large mixing pattern [H] would clearly provide a good leading 
order approximation for the current neutrino mixing pattern, if the atmospheric neutrino 
mixing angle 6*23 turns out to lie in the first octant. However, the second octant of 6*23 can 
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not be ruled out and is supported by the analyses in Refs. [H1IS]- In this case the texture 

sin #12 ~ A, sin# 13 ~ A, sin #23 ~ 1 — A (2) 

could be taken as a viable model building standard. We shall refer to two mixing patterns 
as BLi and BL 2 textures respectively. The difference between BLi and BL 2 mixing lies in the 
order of magnitude of the atmospheric mixing angle 623, the BLi mixing pattern would be 
favored if future experiments establish that 823 belongs to the first octant and the deviation 
from maximal mixing is somewhat large, otherwise BL 2 mixing is preferred. It is well-known 
that the observed hierarchies of masses and flavor mixing in the quarks and charged leptons 
sectors can be conveniently characterized by the Cabibbo angle. As a result the BL t and BL 2 
parameterization may have deep implications for the theoretical formulation of the ultimate 
unified theory of flavor. A lot of work in the literature has demonstrated that the smallness 
and hierarchy of the quark masses and mixing angles can be naturally generated in theories 
which, at low energy, are described effectively by an Abelian horizontal symmetry, which 
is explicitly broken by a small parameter [T7HT5] . It certainly follows a natural path to 
try and apply these ideas on Abelian family symmetries developed for the quarks to the 
lepton sector. In this work, we shall investigate whether and how the BLi and BL 2 textures 
can be reproduced naturally from the Abelian horizontal flavor symmetry. For generality we 
assume that the light neutrino masses arise from lepton-number-violating effective Weinberg- 
like operators. 

The paper is organized as follows. In section |2j we present the effective low energy 
theory for the Abelian U(l) flavor symmetry and its extension to U(l) x Z m x Z„. We 
find that, in order to produce the BL t texture without fine tuning, the family symmetry 
should be U(l) x Z m x Z n with m and n of opposite parity. Models for the BLi an d BL 2 
schemes are constructed in section Eland section H] respectively; these models are extended to 
include quarks within the SU (5) grand unified theory, the observed patterns of both quark 
and lepton masses and flavor mixings are reproduced, and the general phenomenological 
predictions of the models are discussed. Finally, our conclusions are summarized in section 

El 

2 Theoretical Framework 

Our theoretical framework is defined as follows. For definiteness we consider a low energy 
effective theory with the same particle content as the Supersymmetric Standard Model (SM). 
In addition to supersymmetry and the SM gauge symmetry, we introduce a horizontal U(l) 
symmetry and a SM singlet chiral superfield which is charged under the U(l) family 
symmetry; without loss of generality, we normalize its charge to —1. The effective Yukawa 
couplings of the quarks and leptons are generated from nonrenormalizable superpotential 
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terms of the form 

W = {y u ) l3 Q t U c 3 H u l-\ + (y d ) tj Q t D<H d I -J 

+ {y e ) ij L l E c J H d -j + (y v )ij—LiLjH u H u i-\ (3) 

where -Ff U)( j are Higgs doublets, and Lj are the left-handed quark and lepton doublets 
respectively, £/J, -Dj and £^ are the right-handed up type quark, down type quark and 
charged lepton superfields respectively, i, j are generation indices. The parameter A is the 
cutoff scale of the U (1) symmetry, and F(ip) denotes the U(l) charge of the field ip. Note that 
F(H U ) and F(H d ) do not appear in the exponents, since one can always set the horizontal 
charges of the Higgs doublet H u and H u to zero by redefinition of the U(l) charges. The 
last term of Eq.([n]) is the high dimensional version of the effective lepton- number- violating 
Weinberg operator. 

After the horizontal symmetry is broken by the vacuum expectation value (VEV) of the 
scalar component of the superfield 0, one obtains the quark and lepton mass matrices whose 
elements are suppressed by powers of the small parameter (0)/A, which for simplicity is 
usually assumed to be characterized by the Cabibbo angle, i.e., A = (0)/A, then we have 

(M e )« = (yehX^^Vu, (M v )ij = (y v hX F ^ +F ^ V f (4) 

where v U)d = (H Utd ) is the electroweak scale VEV of the Higgs doublet H u ^. The factors 
{l/ujij, (Vdjij, (Ve)ij and (yu)ij are not constrained by the flavor symmetry and are usually 
assumed to be of order one, the freedom in these factors is used in order to obtain a quan- 
titative description of the fermion masses and flavor mixings. Since the holomorphicity of 
the superpotential forbids nonrenormalizable terms with a negative power of the superfield 
0, one has (M u ) tJ = if F(Qi) + F(U^) < 0. Similarly = if F(Q t ) + F{Df) < 0, 

(Me)ij = if F(Li) + F{Ef) < and [M v ) {j = if F{U) + F(Lj) < 0. 

In our framework, the light neutrino masses are generated by the high dimensional effec- 
tive Weinberg operators shown in the last term of Eq.()3]), consequently the light neutrinos 
are Majorana particles and its mass matrix M u is symmetric with {M u )ij = {M u )ji 0. Fur- 
thermore, if all the horizontal charges are positive, the hierarchial structure of the mass 
matrices shown in Eq.(jH) allows a simple order of magnitude estimate for the various mass 
ratios and mixing angles: 

Hhi. ^ x FiQl) ~ F[Qj)+F{u i ) ~ F{u !> ) , — ~ \ F iQi)- F (Qi)+ F i D i)- F i D j) Vi - ~ A F(Qs)_F{Qj) , 
m Uj ' m dj ' v 

H!i ^ X 2[F(L i )-F(L j )]^ ™* _ x F( Ll )-F(L,)+F(E-)-F(E-) ^ gin ^ • ~ A F{Ll)_F(Lj) (5) 

1 Note that if we introduce three right-handed neutrino superfields Nf to generate light neutrino mass 
via type I see-saw mechanism, the structure of the light neutrino mass matrix M v is independent of the 
Nf charge assignment [23H2I], unless there are holomorphic zeros in neutrino Dirac mass matrix Mu or in 
Majorana mass matrix Mm for the heavy fields N c . 
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where is the light neutrino mass, Vij denotes the element of the quark CKM mixing 
matrix. We note that the sign " ~ " implies that there is a unknown order one coefficient 
in each relation, so that the actual value of the mass ratios and mixing angles may slightly 
depart from the naive "power counting" estimate. Moreover, if some fields carry negative F 
charges, then holomorphy plays an important role and the estimates (JHJ) could be violated 
as well. For the BI^ mixing pattern, both sin6* 12 and sin# 23 are of order A, then we should 
require 

F{L X ) = F(L 2 ) + 1, F(L 2 ) = F(L 3 ) + 1 (6) 

This implies F(Li) = F(Lz) + 2, as a result, we have sin #13 ~ A 2 . Therefore we conclude that 
the BLi mixing pattern can not be naturally produced from a pure U(l) flavor symmetry. 
Turning to the BL 2 mixing pattern given by sin 623 ~ 1, sin #12 ~ A and sin #13 ~ A, one 
should choose 

F(L 2 ) = F(L 3 ) = F{In) ~ 1 (7) 



Then we have the (2i) and (3i) (i = 1,2,3) entries of the charged lepton mass matrix 
are of the same order, hence the diagonalization of the charged lepton mass matrix leads to 
large 2-3 mixing. In addition, we obtain 



My ~ X 2F ^ 



( A ' 

A 



A A \ 
1 1 

1 1 / 



K 
A 



(8) 



Clearly the (2-3) sector of the light neutrino mass matrix has a democratic structure, thus 
large mixing in this 23 sector is naturally obtained. However, barring the presence of special 
cancellations, the masses of the second and the third light neutrinos are typically expected to 
be of the same order in this case. As a result the three light neutrinos are quasi-degenerate 
and strong parameter fine tuning is required in order to account for the hierarchy between 
the measured mass squared differences Am 2 ol and Am 2 tm . 



In order to avoid this kind of fine tuning in obtaining and acceptable pattern of neutrino 
oscillation parameters, we must go beyond the pure U(l) flavor symmetry case considered 
above. Let us now move to the extended flavor symmetry U(l) x Z m x Z n . This kind of 
Abelian symmetry is somewhat complex and not yet fully discussed, as far as we know, 
since most previous work concentrated on U(l) or U(l) x Z m flavor symmetry. We now 
consider three SM singlet superfields 61, 2 and 3 with horizontal charges 

6! : (-1,0,0), 2 : (0,-1,0), 6 3 : (0,0,-1) (9) 

which in general can have different vacuum expectation values (©1), (0 2 ) and (©3). For 
simplicity we take in what follows, (0i)/A ~ A, (0 2 }/A ~ A and (Q^/A ~ A. The effective 



i 



Yukawa couplings are given by extending Eq.Q with new flavons 0i, 02 and G3 as follows: 



W = (y u ) ij Q i U*H u ^ 

O \F(Qi)+F(Dj) ,fs \ [Zm(Qi)+Z m (D?)] , fi ^[Zn(Qi)+Z n (D?j\ 



(Ue)ijLiEjH d 



q \ F{Li)+F{E9) /Q\ [Z m (Li)+Z m (E])] / q \[Zn(Li)+Z n (Ef)] 
F{Li)+F{Lj) . 

)] /Q \ [^n(i l )+Z„(L j )] 



+(y l/ ) ijJ L l L j H u H u ^j i^-fj (^j (10) 

where Z m n (^) is the Z m n charge of the field and the brackets [. . .] around the exponents 
denotes that we are modding out by m (n) according to the Z m (Z n ) addition rule, namely: 

r / ,0 r t if r < m , , 

[ZjQi) + Z m {Uf)] = { (11) 

r — m it r > m 

where r = Z m (Qi) + Z m {Uf). We note that the charge assignments of the Higgs doublets if u 
and have been set to be (0, 0, 0) by redefining the flavor symmetry charges of the fields. 
Thus the fermion mass matrix can be expressed in term of the horizontal charges as 

(M u )ij = (^) u .A F ^) +F ^) + [^^) +z ™( f/ P] + [ z "^) +z "(^)]^, 

(Me)* = ( ye ) u .A F ^)+ F ^)+[^^)+^^)] + [^^)+^^)]^, 

(M u )ij = {y v ) . i A F ( i *)+i ;i (^)+[^(^)+^(^)]+[^(i < )+^(^)] ^ . ( 12 ) 
Consider the quark sector, the flavor mixing angles there are given by 

Vf- rsj \( F (Q*)+ F ( DC j))-( F (Qj) + F ( DC 3))+[ Z rn(^ (14) 

For m = n = the CKM matrix elements describing the charged current weak interaction 
of quarks behave approximatively as V^ d ~ \ F (Qi)~ F (Qj) anc l therefore the CKM mixing 
Vckm — ' V d is expected to scale as Vckm^ ~ \ F (Qi)- F (Qj) . In order to compare with 
the pure U{1) horizontal symmetry case, we can define an effective flavor charge in the 
general case m ^ n ^ as, 

F eff (if>) = F(iP) + ZM) + Z n {$) . (15) 

Then it is clear that 

where a, /3 = 0,1 and we have used Eq. (JTTJ and the fact that 

[Z m (Qi) + Z m {Uf)\ - [Z m (Qj) + Z m {Uf)\ = Z m (Qi) - Z m (Qj) ± am (17) 
[Z m (Qi) + Z m {Df)} - [Z m {Q 3 ) + Z m (D])] = Z m (Qi) - Z m {Q 3 ) ± am (18) 
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where a = 0, 1. The condition for the value ±/?n follows similarly. Similarly for the lepton 
sector one obtains 

y}. r^j \Feff(Li)-F ef} iL j )±am±t3n _ 

Therefore the masses and mixing angles can be enhanced or suppressed by \ ±m±n rela- 
tive to the scaling predictions obtained when the family symmetry is the continuous flavor 
symmetry U(l)i x f/(l) 2 x U(l)s because of the discrete nature of Z m x Z n . Note that in 
the case where the light neutrino masses are generated by type I see-saw mechanism and 
all fermion charges are positive, the neutrino masses and mixing angles still do not depend 
on the details of the right-handed neutrino sector, except for the possible enhancement or 
suppression associated to the Z m x Z n flavor symmetry. 



Furthermore, when the flavor symmetry is reduced to U(l) x Z m by taking n = 0, all 
the above results remain valid. It is remarkable that we can employ the U(l) x Z m flavor 
symmetry to maintain the BL 2 mixing while achieving very different neutrino masses without 
fine tuning. We shall restrict our attention to the case of a Z 2 symmetry which is the minimal 
nontrivial Z m group (see for example the explicit model construction given in sec. H] below). 
In this case just the Z m symmetry can reproduce a hierarchy in neutrino masses of order A 2 
consistent with the observed ratio of solar-to-atmopsheric splittings. 

In contrast, note that since the reactor neutrino mixing is necessarily of order sin 6*13 ~ 
\ 2±am the U(l) x Z m flavor symmetry can not produce the BLi mixing pattern. Indeed for 
such BLi texture one has sin# 12 ~ A and sin #23 ~ A, which is in conflict with the required 
linear behaviour of the reactor mixing angle sin #13 ~ A. Note parenthetically that the Z\ 
group consists of only the identity element, so the group U(l) x Z\ is isomorphic to U(l), 
and the Z\ charge of field is 0, hence the flavor symmetry U(l) x Z\ produces a wrong scaling 
behaviour sin $13 ~ A . 

We now turn to the realistic case of the U(l) x Z m x Z n flavor symmetry. If both solar 
and atmospheric neutrino mixing angles are of order A then the reactor angle would be 
constrained to be or order sin 613 ~ ^2±am±/3n_ ^ g a resu }t ; one can have sin #13 ~ A if 
the parity of m and n is opposite. This is an interesting observation of the present work. 
In section [31 a concrete model for the BLi mixing pattern is presented, based on the flavor 
symmetry U(l) x Z 3 x Z 4 . 

Since an Abelian flavor symmetry can not predict the exact value of the 0(1) coefficients 
in front of each invariant operator, we must content ourselves with explaining the orders of 
magnitude of fermion masses and flavor mixing parameters. To identify the phenomenolog- 
ically acceptable mass matrices, we will estimate the various mass ratios and mixing angles 
as approximate powers of the small parameter A. The hierarchies in the quark mixing an- 
gles are clearly displayed in Wolfenstein's truncated form (20] of the parameterization of the 



6 



/ 



v ( 



CKM 



(20) 



CKM matrix p 

1 - A 2 /2 A AX 3 {p-t V ) \ 

-A 1 - A 2 /2 AX 2 

\ AX 3 (1- p-ir]) -AX 2 1 

where the quantities A, p and rj are experimentally determined to be of order one. Therefore 
the order of magnitude of the three mixing angles is given in terms of the A as 

|K,|~A, \V cb \~X 2 , |K fe |~A 3 -A 4 (21) 

The charged fermion mass ratios at the grand unified theory (GUT) scale should satisfy [22] 



A J 



as well as 



m v 
m c 

~ A 

m s 

m e 2 
~ A 



A 3 - A 4 , 
A 2 , 



m t 

~ A 



m b 3 
~ A 



m t 



(22) 



(23) 



for the intrafamily hierarchy. The first identity is the well-known b — t unification relation. 
For the neutrinos, we required that the lepton mixing is of BLi or BL 2 type depending on 
the octant of #23 . We note that the exact power of the A for each of the mass ratios and 
mixing angles may vary a little, depending on the value of tan/3 = v u /vd and on the exact 
value one chooses for A. Moreover, the above order of magnitude estimates depend on the 
scale, since all the parameters would run under the renormalization group evolution. In this 
work, all the explicit models are properly constructed to meet the requirement m t /v u ~ 1 
and m b /vd ~ A 3 so that tan/3 ~ 1. As a result, the scale dependence of these mass ratios 
and mixing angles is weak and hence will be neglected in our framework. 



3 Model for BLi mixing 



As has been shown in the previous section, one can reproduce the BLi texture within the 
framework of U(l ) x Z m x Z n family symmetry, where m and n should have different parity. 
For concreteness, we shall use m = 3 and n = 4 for our model. For such symmetry choice 
the possible model realization of the BLi texture is not unique. As a concrete example, here 
the horizontal charges of the lepton fields are taken to be 



Li: (4,1,3), L 2 : (3,2,2), 
E{: (3,2,2), E c 2 : (1,2,2), 

One immediately obtains the charged lepton mass matrix 

/ A 8 A 6 A 8 \ 



El 



(1,1,1), 
(0,0,0). 



M, 



A 7 A 5 A 7 
7 A s ^3 



(24) 



(25) 
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which yields the mass ratios 



I^A 3 , ^~A 2 , (26) 



that are consistent with the experimental requirements. For the charged assignments in 
Eq. (I24p . the light neutrino mass matrix is given by 



My ~ 



/ A 12 A 8 A 7 \ 

A 8 A 7 A 7 
V A 7 A 7 A 6 / 

It predicts the light neutrino mass eigenvalues as follows 



v 2 



A ' < 27) 



mx-A 8 ^, m 2 ~A 7 ^, m 3 ~ A 6 ^ (28) 

The neutrino mass spectrum is normal hierarchy. It is remarkable that this model gives rise to 
m 2 /m 3 ~ A and Am 2 ol / Am 2 tm ~ A 2 , which is in excellent agreement with the experimental 
data. In conventional U(l) or U(l) x Z m flavor symmetries, if any ratio between neutrino 
masses is an odd power of the small breaking parameter, generally the mixing angle between 
the two neutrinos will vanish [23]. The crucial point is that the element (M u )22, which would 
have been 0(A 14 ) under the continuous ^(1)1 x £^(1)2 x ^(1)3 symmetry, is enhanced to 
0(X 7 ) due to the discrete symmetry Z 3 x Z 4 . Diagonalizing the mass matrices in Eq. (l25|) 
and Eq.( 127|) by the standard perturbative techniques described in Refs. [IS1I231I2S], we get 
the three lepton flavor mixing angles 

sin# 12 ~A, sin# 13 ~ A, sin6 l 2 3 ~A. (29) 

Hence the BLi pattern is produced automatically. Note that the solar neutrino mixing 
sin 612 arises from order A contributions from the diagonalization of both M e and M u , while 
at leading order the reactor and the atmospheric neutrino mixing angles receive contribution 
only from the neutrino mass matrix M u . The off-diagonal elements (M I/ ) 13 and (M I/ ) 23 are 
enhanced by Z4 and Z 3 respectively, hence we have sin#i 3 ~ A and sin #23 ~ A instead of 
the naive expectations sin 813 ~ A 5 and sin 6*2 3 ~ A 4 characteristic of the continuous flavor 
symmetry case. 

In the following, we shall extend the model to encompass also quark sector. Since grand 
unified theories (GUTS) relate quarks and leptons, the transformation properties of quark 
fields can be determined from those of leptons. In order to give a successful description of 
the observed fermion mass hierarchies and mixings simultaneously under the same flavor 
symmetry acting on quarks and leptons we work in the framework of SU(5), for definiteness. 
Another motivation of considering SU(5) unification is the anomaly cancellation. If the 
U(l) flavor symmetry is gauged then a general assignment of flavor charges to the fields 
will be anomalous. One can imagine the anomaly to be canceled via the Green-Schwarz 
mechanism [26J, however one must check whether the correct relations are satisfied [27]. A 
convenient way to ensure that the flavor charges are amenable to cancellation is to have the 



flavor symmetry to commute with the SU(5) group [2B]. In the conventional SU(5) grand 
unified theory, the fields D\ and Lj of the same generation are assigned to a 5 multiplet, 
the fields Qi, U? and Ef are unified in the 10 representation. Since the flavor symmetry is 
required to commute with the gauge symmetry, this means that the fields in the each gauge 
multiplet transform in the same way under the flavor symmetry. Consequently the quantum 
numbers of the quark fields under the flavor symmetry U(l) x Z 3 x Z 4 are as follows: 

Q L1 : (3,2,2), Q L2 : (1,2,2), Q L3 : (0,0,0), 

U{ : (3,2,2), U c 2 : (1,2,2), U c 3 : (0,0,0), 

D\: (4,1,3), D c 2 : (3,2,2), D c z : (1,1,1). (30) 

We note that although there are many possible assignments to produce the BLi texture 
in the neutrino sector, only a few of them can satisfy the quark sector phenomenological 
constraints within SU(5). It is well-known that the minimal SU(5) grand unified theory 
predicts that the down type quark mass matrix is the transpose of the charged lepton mass 
matrix, therefore the down type quarks and charged lepton masses are closely related : 
n^e = n^d, — m s and m T = mb, which are in gross disagreement with the measured 
fermion masses and must be corrected [29]. This can be done through the contribution of 
renormalizable [2H] or non-renormalizable [30] operators to the Yukawa matrices. Following 
Ref. [HT], we introduce an additional U(l) x Z 3 x Z 4 singlet superfield E transforming as a 
75 of SU(5), which has non-renormalizable couplings to fermions of the form 5 lOiJg E/A. 
The Yukawa couplings of the down type quark and charged leptons then arise from the two 
SU(5) x U(l) x Z 3 x Z 4 invariant superpotential terms q 

T 

HtJ • (31) 

which, after the scalar components of E acquires a VEV, lead to: 

(Y^ = {{C^ + KiC^X^+^M^^ , 

(Y e )ij = ((CJij - 3k (c-^.^A^CQ^+^^^CQ^+^^l+^CQ^+^^f)] ? (32) 

where k = (E)/A, which breaks the transposition relation between Y^ and Y e , and can 
explain the difference between down type quark and charged lepton masses. In our numerical 
fits, we take k = 0.3 for illustration, and find that realistic values for down-type quark and 
charged lepton masses can be reproduced. The superpotential for the up type quark mass is 



W d = (iQ^CdifijHi + 1 10,(^)^5^5) (jf) 



Q X \Z A {Wi)+Z A {5j 



Q\F(10i)+F(10j) /Q\ [Z 3 (10i)+Z4(10jj\ /@ \ [2 4 (10 i )+Z 4 (10 



W u = 10 i (C a ) id 10 j H s [ 1 ^\ f-^l (^) , (33) 



2 The 75 could in principle also give a contribution in the up sector namely. However, following Ref. pi] 
we neglect such a term since it is not needed to reproduce the up-type quark masses. 
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where one has (C 3 )jj = (C 3 )jj due to the constraint of the SU(5) gauge symmetry. Then 
one can express the effective Yukawa couplings for the up type quark in terms of the flavor 
symmetry charges as 



(y„: 



ij 



^Q. i y..)F{Q l )+F{Q J )+[z. i {Q l )+Z 3 {Q ] )] + [Z i (Q l )+Z4Q ] )] 



(34) 



With the assignments dictated by Eq.f pOj) . one has the following patterns for the up and 
down type quark mass matrices, 





( A? 


A 5 


A 7 








( A§ 


A 7 


A 7 


\ 




A 5 


A 3 


A 5 




v u , 


M d ~ 


A 6 


A 5 


A 5 






U 7 


A 5 


1 


) 






U 8 


A 7 


A 3 


/ 



(35) 



which yields 



IK. 



A^ 



cb\ 



ub 



A . 



(36) 



We note that both the up and down quark sector contribute A 2 to the mixing element 
\V US \, therefore an accidental enhancement of 0(X~ l ) among the undetermined order one 
coefficients (Ci)y, (C^y and (Cs)ij is required in order to describe the correct Cabibbo 
angle. The remaining CKM mixing angles \V c b\ and \V u b\ arise solely from the diagonalization 
of the down type quark mass matrix M d . In addition, the pattern given by Eq. (l3"o"|) leads to 
the following quark mass scalings 



m u ~ X 7 v u , m c ~ X 3 v u , m t ~ v u , 
m d ~ X 8 v d , m s ~ X 5 v d , m b ~ X 3 v d . 



(37) 



which describe the experimental data satisfactorily. Note that the second term in Eq. fl3ip 
accounts for the mass difference between the down-type quarks and charged leptons, allowing 
for an acceptable charged fermion mass pattern. 

In order to see in a quantitative way how well the model describes the observed values of 
the fermion masses and mixings, we perform a numerical analysis, within three independent 
different seeding methods, namely flat, gaussian and exponential distributions. Namely, the 
modulus of the undetermined order one coefficients are taken to be random numbers with fiat, 
gaussian and exponential distributions in turn, the corresponding phases are varied between 
and 2tt. The probability density function f(x) of the three distributions is well-known 



/(*) 



l 

b—a 



a < x < b 
x < a ot x > b 



(38) 



for flat distribution, we take a 
case of gaussian distribution, 



1/3 and b = 3 for illustration in the present work. In the 
f{x) = ■ (39) 



2na 
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We set the mean [i = 1 and the standard deviation cr = 1.5 in our numerical calculation. 
The probability density function for the exponential distribution is 



Its statistic mean is 1/A, and A is taken to be 1 as a typical value for numerical simulation. 
To the extent that our results are independent of the choice of seeding method, they are 
robust and not simply an artifact of the choice of the seed function. 

The coefficients (Ci)jj, (C 2 )ij, (C 3 )ij and (y u )ij are treated as random complex numbers 
with arbitrary phases and absolute value in the interval of [1/3,3]. Then we calculate the 
quark and lepton masses as well as the CKM and lepton mixing matrix entries which are 
required to lie in the experimentally allowed ranges. The numerical results are found to 
be nicely consistent with the above theoretical estimates and qualitative discussions. Since 
the flavor parameters of the quark sector are precisely measured, here we focus on the 
neutrino sector. As an example the predicted distributions for the light neutrino masses and 
atmospheric mixing parameter are shown in Fig. [TJ The the light neutrino masses follow 
the normal hierarchy pattern and, for all the points produced, though all non-vanishing, 
they are rather tiny, with most of the expected mi values below 0.015 eV. As to the 
mixing angles, no specific values of 6 , 12 and #13 are favored within 3a, and hence they are 
not shown in the figure In contrast, however, the atmospheric neutrino mixing angle #23 
obeys sin 2 #23 < 1/2, which means that non- maximal #23 values are preferred, as indicated 
by current neutrino oscillation global analyses post-Neutrino 2012 [6HH], with a preference 
for the first octant. This has been one of our motivations for introducing BL t mixing pattern, 
which leads to sin 6*23 of order A at leading order. 

The rare process neutrinoless double beta decay (0z/2/3) constitutes an important probe 
for the Majorana nature of neutrino and lepton number violation [32], a sizable number of 
new experiments is currently running, under construction, or in the planing phase. The 
histogram for the distribution of the effective 0z/2/3-decay mass \m ee \ and its correlation with 
the lightest neutrino mass mi are given in Fig. [2j We also show the future sensitivity on 
the lightest neutrino mass of 0.2 eV from the KATRIN experiment [33]. The horizontal lines 
represent the sensitivities of the future Oz/2/3-decay experiments CUORE [M] and MAJO- 
RANA [35j/GERDA III [35], which are approximately 18 meV and 12 meV respectively. 
Clealry the expected effective mass \m ee \ is predicted to be far below the sensitivities of 
the planned 0u2(3 experiments. The reason for this is the strong destructive intereference 
amongst the three light neutrinos, as seen in the right panel. As a result, if 0z/2/3 decay will 
be detected in the near future, our construction would be ruled out. 

3 Similarly, we can hardly see any specific preferred pattern for the CP violating phases <5, ip\ and c/?2, 
hence, as before, these are not shown. 




(40) 
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Figure 1: Histograms for the distribution of light neutrino masses and atmospheric neutrino mixing parameter 
in the BLi model. In the second row, left, middle and right panels are obtained with using different seed 
procedures for the order one Yukawa coefficients, namely flat, exponential and Gaussian, respectively, from 
left to right. 
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Figure 2: Histogram of the effective mass |m ee | (left panel) and the scatter plot of \m ee \ versus the lightest 
neutrino mass mi (right panel) for the BLi model. The colored bands represent the regions for the 3<r ranges 
of the oscillation parameters in the normal and inverted neutrino mass spectrum respectively. The future 
sensitivity of 0.2 eV of KATRIN experiment is shown by the vertical solid line, while the future expected 
bounds on \m ee \ from CUORE and MAJORANA/GERDA III experiments are represented by horizontal 
lines. 



4 Model for BL 2 mixing 

As explained in section |2l the order one atmospheric neutrino mixing sin #23 ~ 1 generi- 
cally implies that the corresponding masses of ^2 and z/3 are of the same order of magnitude 
within pure U(l) family symmetry schemes. As a result the neutrino mass spectrum is 
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quasi-degenerate and strong fine tuning is required in order to account for the measured 
mass-squared differences Am^ and Am 2 tm . Furthermore, the renormalization group evo- 
lution effects could drastically enhance the neutrino mixing angles due to the degeneracy, 
so that the BL 2 texture would be spoiled at the electroweak scale. This can be avoided by 
extending the flavor symmetry to U(l) x Z m . Now the whole flavor symmetry is chosen to 
be U(l) x Z 2 , the lepton fields carry the following U(l) x Z 2 charges: 

L 1 : (3,0), L 2 : (3,1), L 3 : (2,0), 

E{: (4,0), E c 2 : (2,1), E c 3 : (0,1), (41) 
Then the light neutrino mass matrix is given, apart from the order one coefficients, as 

/ A 6 A 7 A 5 \ 
M v ~ A 7 A 6 A 6 
V A 5 A 6 A 4 / 

which yields 

mi~A 6 ^/A, m 2 ~A 6 ^/A, m 3 ~ AS 2 /A (42) 

One sees that the first two light neutrinos are quasi-degenerate in this model, and their masses 
are suppressed by C(A 2 ) with respect to the third one ggg. This prediction is consistent 
with the observation that the solar neutrino mass difference Am^ is much smaller than the 
atmospheric neutrino mass difference Am^ tm . On the other hand, the charged lepton mass 
matrix takes the following form 



A 



/ A 7 A 6 A 4 \ 

A 8 A 5 A 3 
V A 6 A 5 A 3 / 



v d , (43) 



which has a "lopsided" structure, a large 2-3 mixing arises from the diagonalization of M e . 
Obviously it also give the correct order of magnitude for the charged lepton mass ratios. 
Combing the contribution from both the neutrino and the charged lepton mass matrices 
diagonalization, the leptonic mixing angles are given by 

sin# 12 ~ A, sin# 13 ~ A sin# 2 3 ~ 1 • (44) 

This is exactly the desired BL 2 mixing pattern, eq. fl2]). Here we would like to point that, 
sine the Super-Kamiokande data indicted large the atmospheric neutrino, perhaps even max- 
imal [38], there have been several attempts to account for the large atmospheric neutrino 
mixing sin 6*23 ~ 1 in terms of Abelian flavor symmetries [37]. However, it was usually as- 
sumed that the reactor angle #13 was rather small, at most of order A 2 at that time [39] . In 
contrast in our construction the consistency between large sin 6*23 and sizeable sin # 13 mixing 
angles emerges naturally. 
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In what follows, we extend the model to include quarks within and SU(5) unified frame- 
work. The fields Qi and U? together with Ef within the same generation fill out the 10 
representation, while Df and the left-handed lepton doublet Li make up the 5 representa- 
tion. As a result, we can determine the transformation properties of the quark fields under 
the U(l) x Z2 flavor symmetry as follows 



Qi ■■ 


(4,0), 


Q2 : 


(2, 1), 


Q-3-- 


(0,1), 




(4,0), 


U c 2 : 


(2,1), 




(0,1), 


D\: 


(3,0), 




(3,1), 


m 


(2,0) 



(45) 



The up and down quark mass matrices can be determined in as straightforward way, as 
follows 



(46) 





/ 


A 8 


A 7 


A 5 








( A? 


A 8 


A 6 










A 7 


A 1 


A 2 




V u , 


M d ~ 


A 6 


A 5 


A 5 








V 


A 5 


A 2 


1 


) 






U 4 


A 3 


A 3 


J 





which lead to 



\V US \ < 



A, \V cb \ 



^ 2 , \Vub\ 

A 2 , *± - 

TO JL 



- A 3 , 

y2 mb 



A 3 



(47) 



which are in excellent agreement with observed quark mass hierarchies and CKM mixing an- 
gles. As in section [31 we perform a numerical simulation of the expected neutrino oscillation 
parameters. In Fig. [3] we display the resulting histograms for the neutrino mass eigenval- 
ue^. As expected on the basis of the qualitative estimate in Eq. ( l4"2j) . the light neutrino 
mass spectrum is normal hierarchy, the degenerate spectrum being strongly disfavored, and 
almost all the generated points lie in the region of the lightest neutrino mass mi smaller 
than 0.015 eV. The neutrinoless double beta decay predictions are shown in Fig. HI One sees 
that, in contrast with the BLi case, althouth the effective mass \m ee \ is also quite small, with 
\m ee \ around 5 meV preferred, there is a small portion of the parameter space of the model 
where the predictions for \m ee \ approach the future experimental sensitivities. Therefore, 
upcoming high sensitivity 0^2/3-decay experiments will provide a test of the present model. 



5 Conclusions 

The recent neutrino oscillation experimental highlights: (i) rather large value of reactor 
mixing angle #13 and (ii) indication of significant deviation of the atmospheric neutrino 
mixing angle 823 from maximality may change our theoretical approach for constructing 
neutrino mass models. In this paper, we have studied the Wolfenstein-like mixing schemes: 

4 Insofar as the neutrino mixing angles 9ij and CP phases 6, ipi and (fi2 are concerned, we do not obtain 
any special predicted pattern, hence the results are not displayed. 
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Figure 3: Light neutrino masses in the BL 2 model. 
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Figure 4: Histogram of the effective mass |m ee | (left panel) and the |m ee | versus the lightest neutrino mass 
mi correlation (right panel) predicted in the BL 2 model. 



BLi mixing, in which sin 12 ~ A, sin 13 ~ A, sin 2 3 ~ A, and BL 2 mixing, in which sin 12 ~ A, 
sin ~ A, sin 2 3 ~ 1- The largish 13 can be naturally accommodated in both of them, the 
two mixing patterns differ in the order of magnitude of sin 023, the BLi texture is favored for 
023 in the first octant, while BL 2 is preferred for the second octant 023. In order to produce 
the BLi mixing without invoking unnatural cancellation, the Abelian flavor symmetry should 
be Z7(l) x Z m x Z n with the parity of m and n being opposite. A concrete model based 
on U(l) x Z% x Z4 family symmetry is constructed, where the light neutrino mass hierarchy 
m 2 /m 3 ~ A is realized due to the discrete nature of Z 3 x Z 4 . The ratio Am^ ol / Aml tm is 
expected to be of order A 2 in this model, which is in good agreement with experimental data 
in contrast with conventional U (1) or U (1) x Z m flavor symmetry constructions. Furthermore, 
the model is embedded into the SU (5) grand unified theory to describe the quark masses and 
mixing simultaneously. As for the BL 2 mixing, it can be reproduced within the framework 
of pure U(l) flavor symmetry, however, the light neutrino mass spectrum is expected to be 
quasi degenerate, hence fine-tuning of the neutrino mass parameters is needed in order to 
achieve the observed mass-squared differences. To improve upon this situation, the family 
symmetry is enlarged to U(l) x Z2, which gives rise to both large atmospheric neutrino 
mixing sin 6*23 ~ 1 and hierarchical neutrino masses. The model is extended to SU (5) grand 
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unified theory as well. 

We show that both models can give a successful description of the observed quark and 
lepton masses and mixing angles, and the numerical results are nicely in agreement with the 
theoretical estimates and the qualitative discussions. The light neutrinos are normal mass 
hierarchy in both models, quasi degenerate spectrum is strongly disfavored. The present 
framework can not predict the CP violating phases 5, cpi and y?2- The Oz/2/3-decay effective 
mass \m ee \ is predicted to be rather small in both constructions, substantial part of the 
data are below the sensitivity of future experiments except for a region of the BL 2 model 
indicated in Fig.HJ Therefore future 0^2/3-decay experiments such as CUORE, MAJORANA 
and GERDA III will provide important tests of the present models. 
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